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ABSTRACT 

Context. Flow-induced instabilities are relevant for the storage and dynamics of magnetic fields in stellar convection zones and possibly also 
in other astrophysical contexts. 

Aims. We continue the study started in the first paper of this series by considering the stability properties of longitudinal flows along magnetic 
flux tubes. 

Methods. A linear stability analysis was carried out to determine criteria for the onset of instability in the framework of the approximation of 
thin magnetic flux tubes. 

Results. In the non-dissipative case, we find Kelvin-Helmholtz instability for flow velocities exceeding a critical speed that depends on the 
Alfven speed and on the ratio of the internal and external densities. Inclusion of a friction term proportional to the relative transverse velocity 
leads to a friction-driven instability connected with backward (or negative energy) waves. We discuss the physical nature of this instability. In the 
case of a stratified external medium, the Kelvin-Helmholtz instability and the friction-driven instability can set in for flow speeds significantly 
lower than the Alfven speed. 

Conclusions. Dissipative effects can excite flow-driven instability below the thresholds for the Kelvin-Helmholtz and the undulatory (Parker- 
type) instabilities. This may be important for magnetic flux storage in stellar convection zones and for the stability of astrophysical jets. 



• £h 1. Introduction 

><■ 

■ The stability of magnetic structures with field-aligned flows is 
?3 °f potential astrophysical relevance in relation to magnetic flux 
storage in stellar convection zones, for siphon flows and other 
flows in stellar atmospheres and envelopes, as well as for jets 
and collimated outflows . Whil e the first paper of this series 



necti on with resonant a b sorption of wav e s in sh ear layers, 



by iHollweg et all dl990l). iGoossens et all (I1992I) . iTirrv et al 
( 1998). and lAndries & Goossensl (bOOlalbh . 



dSchiissler & Ferriz-M as, 2007) studied the effects of perpen 



dicular flows on the stability of magnetic flux tubes, here we 
consider the case of parallel flows, directed along the tubd3- 

Linear stability properties of field-aligned flows in ideal 
MHD have been considered previously with regard to 



Chandrasekhatl 1 1 96 ll IParketi 1 19641 kavl 198l|; 


Ferrari et al., 


1981; Satya Narayanan & Somasundaram, 1982 


; Rae, 


1983; 


D'Silva&Choudhuril Il99ll Ichenal. U994: Kolesnikov et all 


2004). Dissipati 


ve effects leading to instability in the 


pres- 



The effect of a longitudinal flow a long a thin toroidal mag- 
netic flux tube has been considered by van Ballegooiienl ( 1 1983b 
and Fe rriz-Mas & Schussler dl993ll 19951) in the case of a rotat- 
ing star and by lAchterberd dl996bl) in the case of a Keplerian 
accretion disk. A general formalism for the stability of station- 
ary flows a long curved th in magnetic flux tubes has been de- 
veloped by lSchmittl dl998l) . 

The work presented here is motivated by the unexpected 
appearance of instabilities in numerical simulations of thin 



ence of longitudinal (shear) flows, connected with the ap- 
pearance of negative energ y waves, hav e been discussed by 
Ryutovi dl988i ll990l) and Ijoarder et all d 1997b and, in con- 



1 For flows that do not vary along the tube, this covers both flows 
inside and outside the flux tube, since the two cases are related to each 
other by a simple Galilean transformation of the frame of reference. 



flux tubes in stellar convection zones dHolzwarthi 120021) in 
para meter domains where linear analysis predicts stability 
(e.g.. lFerriz-Mas & Schusslen , ll995h . The instabilities could be 
traced back to the aerodynamic drag force included in the sim- 
ulations, which does not appear in linear analysis owing to its 
quadratic dependence on the relative velocity between flux tube 
and environment. Such instabilities affect the storage of mag- 
netic flux at the bottom of a stellar convection zone, so that 
a systematic investigation of the underlying physical mecha- 
nisms and the consequences for flux storage is called for. Our 
study is carried out in the framework of the thin flux tube 
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approximation dSpruin, 119811) . which is relevant for magnetic 
structures in the deep layers of a stellar convection zone. This 
approximation is still the only possibility to carry out numeri- 
cal simulations of the dynamics of magnetic structures in these 
layers with realistically high values of the plasma /3 (ratio of 
gas pressure to magnetic pressure), since the correspondingly 
low density, temperature, and pressure contrasts are difficult to 
maintain against numerical diffusion in 2D/3D MHD simula- 
tions. 

Of particular importance with respect to the question of flux 
storage is the inclusion of stratification effects, so that the com- 
bined effects of flow-driven instability and Parker instability 
can be studied. We provide a unified linear treatment of the 
Kelvin-Helmholtz instability, dissipative (friction-induced) in- 
stability, and the undulatory (Parker-type) instability in a grav- 
itationally stratified medium. While we restrict ourselves here 
to straight flux tubes in a plane-parallel stratification in order 
to elucidate the physical mechanisms, the subsequent paper in 
the series will consider the case of toroidal flux tubes and will 
quantitatively discuss the relevance of these instabilities in ro- 
tating stars and the dependence on the stratification parameters, 
with particular attention to the problem of flux storage in stellar 
convection zones. 

The paper is organized as follows. We first consider the sta- 
bility of a straight flux tube with longitudinal flow in a (non- 
stratified) uniform background medium in Sect. [2] The treat- 
ment is extended in Sect. 0to include the effect of gravitational 
stratification. Section|4]gives our conclusions. 

2. Uniform background medium 

As a first example, we consider a straight, untwisted magnetic 
flux tube, which harbors a flow along the field lines with con- 
stant velocity, uq, and is embedded in a uniform and static en- 
vironment. Under these circumstances, the equilibrium for a 
thin flux tube is simply given by the condition of total pressure 
balance, 



PeO 



(1) 



where p e o and are the gas pressures in the environ- 
ment and within the flux tube, respectively, and Bo 
const, is the magnetic field strength. We indicate equi- 
librium q uantities by an index "0". Fo l lowing previous 
work (e.g.. ISpruit & van Ballegooiienj. Il982t ISchiissleil Il990t 

I: IschmittL 



Ferriz-Mas & Schiisslerl Il993l Il995b ISchmitA 1 1998h . 



we con- 
sider a Lagrangian displacement, of the equlibrium flux tube, 
expressed in terms of the Frenet basis, {e t o,e n o,ew} in the tan- 
gential, normal, and binormal directions of the equilibrium flux 
tube, 



b0 ■ 



(2) 



We define the direction of e t o such that uq > 0. For a straight 
flux tube, e„o and et,o can be defined as any pair of orthogo- 
nal vectors that are perpendicular to e t o such that {e t o,e n o,£bo) 
forms a right-handed basis. The displacement components are 
functions of time, t, and arc length along the equilibrium flux 
tube, so- 



2.1. Kelvin-Helmholtz instability 

In the ideal case without friction between the flux tube and 
its surroundings (or other dissipative processes), a uniform, 
straight flux tube can only become unstable owing to the 
Kelvin-Helmholtz instability. The linearized equations for the 
three components of the displacement vector are all decoupled 
from each other. For a transverse dis placement in the normal 
direction we have (see Schmitti ll998l) 



//| n + 2m 4' + ("0 - «a)&' = ' 



(3) 



where ua = Bo/(4npn)) 1 ^ 2 is the Alfven velocity. Dots indicate 
(partial) time derivatives and primes denote derivatives with re- 
spect to arc length, so. The factor p = 1 +(p e o/Pio) takes into ac- 
count the co-acceleration of the surrounding medium as a con- 
sequence of a transverse acceleration of the flux tube, lea ding 



to enhanced inert i a (e.g., iSpruin, 11981b iMoreno-Insertis et al 



1996; Achterberg, 1996a). Here, p e o and p;o are the gas densi- 
ties outside and inside the equilibrium flux tube, respectively. 
The corresponding equation for the binormal component of the 
displacement is analoguous to Eq. (0. The tangential com- 
ponent describes (Doppler-shifted) longitudinal tube waves, 
which are of no further interest in this context. 

Inserting into Eq. (0 an exponential ansatz of the form 
£ n oc exp i{kso - ut), with real wavenumber k > and com- 
plex frequency a>, leads to the dispersion relation 



peo - 2uoku> + k Im — u A j = . 



(4) 



Instability, i.e., an exponentially growing mode, corresponds to 
a positive imaginary part of co. In the case k — 0, Eq. (@]l only 
has the trivial solution, while for k + we have instability if 
the flow velocity surpasses a critical velocity, viz. 



MO > 



P 



1/2 



MA ■ 



(5) 



Note that a finite value for the critical velocity requires that 
p > 1, i.e., that the acceleration of the surrounding medium 
is taken into account. In order to identify the character of 
this instability, we compare the condition derived here with 



the res ult for a 'thick' flux tube obtained by Kolesnikov et al 



(120041) for incompressible perturbations. In the thin-tube limit, 
kR <k 1, where R is the tube radius, the dispersion relation 
[their Eq. (27)] becomes identical to our Eq. ©, if we take into 
account the sign flip in their definition of a>. Consequently, our 
analysis provides the thin-tube limit for the criterion for the 



Kelvin -Helmholtz instability. Furthermore, Kol esnikov et al 



(2004) have shown that the critical velocity in the case of a 
thick tube becomes at most by about 10% higher than the 
value determined from Eq. (0, which therefore can be con- 
sidered as a reasonable estimate (lower limit) for all flux tubes. 



Similar criteria ha ve been given bv Ryutoval (fT990) and also by 
Achterbergl(ll982h and by lOsin et alJdl999h . but the connection 



to the Kelvin-Helmholtz instability was not made in the latter 
two papers. 

We have verified further (see Appendix [A} that the result 
does not depend on the chosen frame of reference, i.e., whether 
one considers a flow within the flux tube with the exterior at 
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rest or a (uniform) external flow along a stat ic flux tube. This 



is obvious for the 'thick' flux tube treated by Kol esnikov et al 



(2004), where the perturbation of the exterior is consistently 
taken into account, in contrast to the thin-tube case, where the 
effect of the exterior is solely described by the enhanced inertia 
term. 

2.2. Friction-induced instability 

Up to here, we have assumed the ideal case of a dissipa- 
tionless plasma. Including dissipative processes can lead to 
qualitatively different behaviour. For magnetic flux tubes in 
stellar convection zones, the drag resulting from a motion of 
the tube perpendicular to its axis relative to the surrounding 
medium is an important dissipative process. For flows with 
large Reynolds number, the drag force is roughly prop ortional 
to the square of the rela tive perpendicular velocity (e.g. JParken 
1 975t iBatchelon, 1 1 9671 §5.1 1). In the case of a finite perpendic- 
ular velocity in the undisturbed flux tube equilibrium, a linear 
stabilit y analysis has been carri ed out in the first paper of this 
series (Schussler & Ferriz-Mas, 2007). In the case of a static 



background medium considered here, the drag force is of sec- 
ond order in the perturbations, so that its effect cannot be stud- 
ied by way of a linear stability analysis. We therefore assume a 
Stokes-type friction, which is linear in the perpendicular veloc- 
ity. It turns out that the resulting stability criteria are indepen- 
dent of the value of the proportionality factor (provided that it 
differs from zero) and we have reason to assume that the crite- 
ria are, in fact, independent of the nature of the dissipation pro- 
cess altogether. For instance, numerical thin-tube simulations 
including a quadratic drag force confirm the criteria derived un- 
der the assumption of a Stokes-type friction (see Holzwarth et 
al. in preparation, paper III of this series). We write the Stokes- 
type friction acceleration in the form 



as = -a (H re i • e„) e n o 

where a > is a constant parameter and 

"rel — U\ — M e 



(6) 



(7) 



is the relative velocity, i.e., the difference between the perturbed 
internal velocity, M;, and the external velocity, M e . Here we have 
M e = and thus obtain the linearized expression 



"rel = MoCtO +€ + 

= (MO + tt + "O^tVtO + (£n + MoO^nO , 

where 



(8) 
(9) 



is the perturbed unit normal vector. The linearized part of as 
provides an additional term of the form a^ n in Eq. (0, leading 
to a modification of the dispersion relation given by Eq. 



2 2 / 2 2 \ 

/j.a> — 2u()kw + k ( m - u A ) + iaco = . 



(ID 



The resulting condition for instability (positive imaginary part 
of a>) is simply given by 



for k + 0. Comparison with Eq. © shows that instability in the 
case with friction sets in for a lower critical velocity than for the 
Kelvin-Helmholtz instability, mkhi = — 1)«a > u a- This 

suggests that the effect of friction leads to a new kind of insta- 
bility, which does not exist in the dissipationless case, although 
the condition for instability given by Eq. (fT2l is independent 
of the actual value of the friction parameter, as long as a does 
not vanish. In the limit a — > 0, the growth rate of the friction- 
induced instability in the velocity range ma < uo < «khi goes to 
zero, while for uq > mkhi we have Kelvin-Helmholtz instabil- 
ity. The criterion Eq. ( fT2l is also independent of the enhanced- 
inertia parameter, /i, which indicates that the instability mech- 
anism is not closely related to the perturbation of the external 
medium, such as in the case of the Kelvin-Helmholtz instabil- 
ity. 

In order to elucidate the nature of the friction-induced in- 
stability, we consider the phase velocity, M p h = Re(a>/k), of the 
unstable mode. It is readily shown that this mode has M p h > 
for Mq > m a , so that the unstable mode is prograde, i.e., it prop- 
agates in the direction of the flow. For flow velocities below the 
instability limit, the direction of propagation of this mode is ret- 
rograde, against the flow. Modes whose direction of propaga- 
tion becomes reve rsed by the flow are sometimes called back 



ward waves (e.g., iNakariakov & Robertsl 1 19951 Ijoarder etal 



1997b . The reversal of the propagation direction at the insta- 
bility threshold means that both real and imaginary parts of 
io change sign when uq - ma goes through zero. We see from 
Eq. (fTTT i that this corresponds also to a sign change of the 
friction-related term, iau), in the dispersion relation: the effect 
of friction changes from wave damping to wave excitation. 

Backward waves can also be nega tive energy waves (e.g., 
Cairn si 1979 : Lashmore-Davies 2005), which owe their name 
to the fact that they can grow in amplitude when the total ki- 
netic energy of the system (energy of the flow plus that of 
the wave) is dec r eased (e.g., by some dissipative process). 
Following ICairnsI (11979b . we call a given mode a negative en- 
ergy wave if the relation 



C = o)^- < 
oa> 



(13) 



holds, where D represents the left-hand side of the dispersion 
relation in the non-dissipative case, here given by Eq. (0J. In 
general, D has to be multiplied by a factor + 1 such that C > 
in the case of a system without flows (uq = 0). In our case, the 
condition for negative energy waves becomes 



to \Ll) - 



uok 



<0. 



(10) Considering the solutions of Eq. 



uok UQk 



2 \l/2 



+ — f 



(14) 



(15) 



"*0 / 



Uq > Ma 



(12) 



we find that a> + (the prograde mode) never satisfies the con- 
dition given by Eq. <TT~4T > while w_ (the retrograde mode) in 
fact represents a negative energy wave provided that Uq > u\. 
Consequently, we find that the following three statements are 
all equivalent in our case: (1) we have friction-induced insta- 
bility; (2) the system supports a backward wave mode; and (3) 
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the system has a negative energy wave. The fact that our unsta- 
ble mode is a negative energy wave suggests that not only the 
exact value of the friction coefficient is irrelevant for the stabil- 
ity criterion but that even the nature of the dissipation process is 
largely irrelevant. In fact, it has been show n that energy loss by 
radiation of MHD waves (Ryutova, 1988), thermal conduction 



(IJoarder et al.lll997l) . or resonant absorption (ITirry et al.lll998l) 



can provide the necessary dissipation for the development of 
the instability. On the other hand, the fact that Eq. ( fT3l l is not 
invariant with respect to Galilei transformations raises ques- 
tions concerning th e general validity of the concept of negativ e 
energy waves (e.g., Walker, 2000; lAndries & GoossensL 12002 ) 
so that we do not pursue this approach any further. 



"rel = ("i 




Fig. 1. Schematic sketch illustrating the physical mechanism of 
the friction-driven flux tube instability. A leftward propagating 
transverse wave along the flux tube (indicated by the curved 
line) is amplified if the speed of the leftward external flow (in 
the frame of reference where the plasma in the equilibrium flux 
tube is at rest) exceeds the phase velocity of the wave. See main 
text for a detailed discussion. 

The physical mechanism driving the friction-induced insta- 
bility can be most easily understood by considering the frame 
of reference in which the plasma in the equilibrium flux tube 
is at rest. We then have an external flow along the equilib- 
rium flux tube with a velocity « e = -Moe t o- Figure Q] schemati- 
cally shows a snapshot of a segment of the perturbed flux tube 
(black curve), together with velocity and acceleration vectors 
relevant for the effect of friction. We assume a transverse wave 
propagating leftward, i.e., in the direction of the external flow. 
The velocity u\ of a mass element of the tube corresponding 
to that wave motion is vertically upward for the tube segment 
and phase shown. The relative velocity between the tube and 
its environment that is relevant for the friction is given by the 
component of the velocity difference, u\ - u e , along the instan- 
taneous unit normal vector of the mass element, n, cf. Eqs. © 
and @. Depending on the speed of the external flow, this pro- 
jection can be parallel or antiparallel to the normal vector. We 
consider two cases: 1) for a sufficiently slow flow (w et i, blue 
vector arrows in Fig.[T]i, the relative velocity (H re i,i) is paral- 



lel to the normal direction, so that the acceleration by friction 
(ai) is antiparallel to the normal and thus decelerates the up- 
ward motion of the mass element, corresponding to a damping 
of the wave; 2) in the case of a sufficiently fast flow (H e ,2, red 
vector arrows), the relative velocity (n re i,2) becomes antiparal- 
lel to the normal direction, so that the corresponding friction 
acceleration (a;) is now parallel to the normal and accelerates 
the upward motion of the mass element, thus amplifying the 
wave. For a wave propagating to the right (reverse direction of 
Mi), we always obtain damping. 

We can put this in quantitative terms by considering the 
relative velocity, 

"rel = («i - «e) • «n = & - "0^n • (16) 

Inserting a wave solution £ n oc sin(&so _ ^Jt), we readily find 

«rel = («ph + Uo)u p ^n , (17) 

where « p h = a>/k is the phase velocity of the wave. Now, for 
w p h > the relative velocity has always the same sign as £ n , 
so that the wave is damped by friction. The same is true for 
«ph < and uq < -u p h, while for uq > -u p h the sign of the 
factor in front of £ n in Eq. (fTTT i reverses, so that the friction 
force excites the wave. Figuratively, one can imagine the wave 
being amplified by a 'tailwind'. 

In the frame of reference with the external medium at rest, 
the condition uq > -M p h transforms into that for the 'backward 
wave', Mph > 0. 

3. Stratified background medium 

We now consider the effect of a gravitational stratification 
of t he external medium, so that buoyancy-driven instabil- 
ity dSpruit & van Ballegooiien , 1982 ; Ferriz-Mas & Schtisslei , 
1993. 1 19951) is included. For a straight, horizontal flux tube, 
mechanical equilibrium then requires neutral buoyancy, i.e., 
PiO = Pe07 in addition to the pressure balance condition given 
by Eq. ([TJ. We define the normal direction of the equilibrium 
flux tube to be along the (constant) gravitational acceleration. 
The linearized evolution of the three components of the dis- 
placement vector is then governed by the following system of 
equations: 

u 2 

£ + 2w £ + («g - - -rNn = .(18) 

u 1 

p£ n + 2 U0 £' n + [U 2 - U 2 A )£ + of n + — ^£ + WMBV^n = ,(19) 

/i& + 2 U0 ? b + (u 2 - u 2 A ) + a& = (20) 

In deriving these equations^ (see also ISchmittt fl 998). we have 
assumed an ideal gas, adiabatic perturbations, and we have 



2 Describing the effect of enhanced inertia for vertical motion by 
simply introducing the constant factor \i into Eq. (1 19b is a du bious pro- 
cedure in the case of a stratified medium (Achterberg, 1996a). Lacking 
a better treatment, we keep this crude description. The resulting stabil- 
ity criteria are independent of \i; only the growth rate of modes dom- 
inated by the Kelvin-Helmholtz instability of normal (vertical) per- 
turbations depends critically on enhanced-inertia effects, so that these 
numbers should be considered with caution. 
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taken the limit [3 = ^npo/B^ » 1, which reflects the condi- 
tions in stellar interiors. This simplifies the expressions without 
changing the character of the various instabilities. In Eqs. ( fT~8b - 
( l20l i. y denotes the ratio of the specific heats, H p is the pressure 
scale height, a the friction pa rameter, and wmbv is the mag netic 



Brunt-Vaisala frequency (cf. lMoreno-Insertis et al.LI 19920 with 



2Hj 



2 - 






7 







P8 



(21) 



where <5 = V - V a d is the superadiabaticity of the external strat- 
ification. For k — a — and w MBV > 0, the perturbed flux tube 
would perform stable oscillations in the normal direction with 
to - wmbv- For k — 0, w MBV > 0, and a + 0, the tube behaves 
like a damped harmonic oscillator. The binormal component, 
Eq. d20l >. is decoupled from the rest of the system and formally 
identical to the non-stratified case treated in the preceding sec- 
tion, so that the same results apply. 

For the coupled normal and tangential components we 
again use an exponential ansatz of the form £ t; n(so, t) = 
|t, n exp(zfao _ Kt>t) and normalize all frequencies by ku A , writ- 
ing to - a>/(ku A ), wmbv = wmbv/(^a), and a = a/(ku A ). 
Inserting this into Eqs. ( TT~8T > and ( fT9l and keeping the same 
symbols for the non-dimensionalized quantities, we obtain the 
linear system 

[to 2 - 2M A to + [m\ - l)] It + iq$ B = (22) 



[fito 2 - 2M A to + (m a - l) • 



~2 



,]#n-fcft = 0, (23) 



where q = l/(yH p k) is a non-dimensional measure of the 
wavelength and M A = uq/u a is the Alfvenic Mach number. 
This leads to a dispersion relation of the form 



jlOJ A + C3& 3 + C2<b 2 + C\(l) + Cq = 



with 



(24) 



C3 = -2(jj + l)M A + i& , 

c 2 = (ju + 1) (m\ - l) + 4M 2 A - w MBV - 2iM A a , 
ci = 2M A o> MBV - {M\ - l) (4M A - /a) , 

K-OK-i-^mbv)-? 2 - 



CO 



3.1. Parker instability Qi = 1, a = 0) 

In the case yu = 1 and 5 = 0, the trans formation u> — > w - Ma 
leads to a simple biquadratic equation (Schmitt, 19981) . The so- 
lutions show that the flow does not affect the stability proper- 
ties in this case. The criterion for instability (positive imaginary 
part of Co) is given by 



,~,2 



< q 2 - 1 . 



(29) 



which corresponds to the result of ISpruit & van Ballegooijen 
d 1982D . The condition marks the onset of the buoyancy-driven 
undular instability or Parker instability. In fact, the stability cri- 
terion for the Parker instability remains unchanged for fi > 1 . 
Note that, owing to our chosen non-dimensionalisation, the 
wavelength-dependence of the Parker instability is hidden in 
wmbv = w M bv/(^a) and q = l/(yH p k). 




0.5 1.0 1.5 

Alfvenic Mach Number 



Fig. 2. Kelvin-Helmholtz instability in the case with external 
stratification (/i = 2, q = 1, o> MBV = 0.2) without friction 
(a = 0). The full and dashed lines give the real and imaginary 
parts, respectively, of the (normalized) frequency, to, for the 
four solutions of the dispersion relation, Eq. d24l >. as functions 
of the Alfvenic Mach number of the longitudinal flow, M A . 
Instability sets in for M A 0.55, when the real parts of the two 
(mainly) transverse modes merge and the imaginary part of to 
becomes positive for one of the modes. The two modes whose 
real parts cross the zero line at M A 0.31 and M A 1.45, re- 
spectively, represent backward modes. They become unstable 
if friction is included (a > 0). 



2.0 





CD 


1.5 




£ 




(25) 


z 




(26) 


_c 




1.0 


(27) 


(J 
c 

CD 




(28) 


> 
< 


0.5 




0.0 

0.0 



Fig. 3. Critical Alfvenic Mach numbers for the onset of 
the friction-induced instability as a function of the non- 
dimensionalized wavelength, q, for w^ BV = 0.2. For one of 
the two modes that can become unstable backward waves, 
the critical Mach number goes to zero as the wavelength ap- 
proaches the value for the onset of the Parker instability, q = 
(l+w MBV ) 1 '' 2 (indicated by the dashed vertical line). The dotted 
line indicates the case q — 1 shown in Fig. [2] 
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Alfvenic Mach Number 
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3 0.00 
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0.5 1.0 1.5 

Alfvenic Mach Number 



2.0 



Fig. 4. Kelvin-Helmholtz and friction-induced instability for 
the same case as that shown in Fig. [2] now with a = 0.2. The 
real and imaginary parts of the normalized frequency, &>, are 
given by the full and dashed lines, respectively. The upper panel 
shows all four modes and the full range of frequencies, while 
the lower panel gives a detail of the zero crossing of the two 
backward modes. 



3.2. Kelvin-Helmholtz instability (jj. > 1, a = 0) 

For n > 1, we have the possibility of a Kelvin-Helmholtz insta- 
bility (cf. Sect. 12. It , now modified by the presence of stratifica- 
tion. If the flux tube is stable with respect to the Parker instabil- 
ity, i.e., if the instability condition given by Eq. ( f29b is not met, 
Kelvin-Helmholtz instability sets in for a flow exceeding a crit- 
ical Alfvenic Mach number, Ma,c, which varies between zero 
for neutral Parker instability (i.e., w^ BV = q 2 - 1) and infinity 
for w^ BV — > oo. Hence, a flux tube equilibrium approaching 
the stability limit of the Parker instability may become Kelvin- 
Helmholtz unstable for arbitrarily low Alfvenic Mach number. 

The physical reason for this behavior is the fact that near 
the stability limit for the Parker instability, most of the sta- 
bilizing effect of the magnetic tension force has already been 
compensated for by the stratification effects, so that a low flow 
velocity is sufficient to destabilize the system. As an example, 



Fig. |2] shows the real and imaginary parts of the complex fre- 
quency, a>, as a function of M\ for = 2, q - \, a - Q, 
and <I^ BV = 0.2. In the limit H p — > oo (q — > 0, wmbv — * 0, 
non-stratified case), the two mode pairs rep resent the transverse 
and longitudinal tube modes dSpruiill98lh . With stratification, 
longitudinal and transverse perturbations are coupled and the 
modes are of a mixed type. Kelvin-Helmholtz instability sets in 
at Ma — 0.55, when the real parts of the frequencies of one pair 
(corresponding to the transverse mode in the unstratified case) 
merge and an unstable mode (positive imaginary part) appears. 
In the non-stratified case, the critical Alfvenic Mach number is 
significantly larger, namely equal to V2 [cf. Eq. Q]. 



2.0 3.3. Friction-induced instability (a + 0) 



In the case a > 0, the effect of friction adds another possibil- 
ity for instability. Guided by the results discussed in Sect. 12.21 
we conjecture that the instability is connected with backward 
waves, leading to a reversal of the sign of the friction term in 
Eq. (f23t . This conjecture is verified by direct numerical so- 
lution of the dispersion relation, Eq. ( l24t . Fig. [2] shows that 
the real part of the frequency for two modes goes from nega- 
tive to positive sign, indicating the presence of two backward 
waves for suffienctly high values of Ma- Writing a> = (a> r , aii), 
the conditions for the presence of a backward wave becoming 
marginally unstable for a> = (0, 0) are 



dcby 



dM A 
and 
da>. 



> 



w=(0,0) 



> 0. 



(30) 



(31) 



&=(0,0) 



The dispersion relation, Eq.fSH shows that u> = (0, 0) is equiv- 
alent to co = 0, which entails the relation 



^MBV (l - M l) 



(l-M A ) 2 



(32) 



Taking the derivative of Eq. d24l ) with respect to Ma and setting 
a) = (0, 0) gives 



dib 



dM A 



dc 
dM A 



0. 



(33) 



S=(0,0) 

Separating this equation into its real and imaginary parts leads 
to 



du) r 



dM A 
and 

d(b{ 



dM A 



&=(0,0) 



&=(0,0) 



dM A c 2 +c 2 J 



CIS 
CI i 



dw r 



dM A 



6>=(0,0) 



(34) 



(35) 



where ci = (ci r , ci i) and co is real. Inserting the expressions 
for these coefficients given in Eqs. ((27) and (|28] i into Eq. d34l 
leads to a quadratic expression, which is positive as long as 
Ma + 0. Similarly, by using Eq. ((27) and the relation given by 
Eq. (l32l . one shows that -C\ijc\ f > 0, so that both conditions 
for unstable backward waves given by Eqs. (f30l > and (|3TT i are 
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fulfilled. This extends our previous result that friction-induced 
instability occurs in the form of unstable backward waves to 
the case of a stratified medium. 

The two mode pairs described by the dispersion relation 
Eq. (l24l result from the combined action of the restoring forces 
of buoyancy (magnetic Brunt- Vaisala mode) and magnetic ten- 
sion (transverse tube mode) in the system. In both cases, the 
retrograde member of the mode pair can become a backward 
wave, so that we expect two modes to become friction-unstable. 

In fact, as already indicated in Fig. [2] each of the mode 
pairs resulting from the dispersion relation Eq. ( l24b exhibits 
a backward wave, becoming unstable at the critical Alfvenic 
Mach numbers given by the two roots of Eq. > 



Mi 



1 + 



,~.2 
MBV 



1 ± 


,~.4 


1/2- 


(36) 




v '"mbv > 







Note that, as long as they do not vanish, neither the value of 
a nor that of the enhanced interia factor, fi, affect the stability 
criteria. It is only the growth rate of the unstable modes that 
depends on the value of a. 

Figure[3]shows the two branches of critical Alfvenic Mach 
numbers resulting from Eq. d36b for w^ BV = 0.2 as a func- 
tion of the non-dimensionalized wavelength, q. For q — > 0, 
i.e. wavelength much smaller than the scale height, the lower 
branch corresponds to the result for the non-stratified case 
(transverse mode) with Ma - — > 1, while the upper branch orig- 
inates from the longitudinal tube mode. For finite wavelength, 
the upper branch goes to larger critical Ma, while the values 
decrease for growing q in the lower branch. At the threshold 
for the onset of the Parker instability, q = (1 + oiL^y) 1 (indi- 
cated by the dashed vertical line in Fig. [3]), Ma,- — > 0. Similar 
to the case of the Kelvin-Helmholtz instability discussed in 
the previous subsection, a growing part of the stabilizing ef- 
fect of the curvature force is compensated by buoyancy effects 
as the wavelength approaches the critical wavelength for the 
Parker instability, so that lower flow speeds suffice to trigger the 
friction-induced instability. The dotted vertical line in Fig.[3]in- 
dicates the case q — 1 shown in Fig. [2] the intersections of this 
line with the two branches give the Mach numbers at which the 
corresponding two modes become backward waves (i.e., where 
the real part of a> changes its sign). 

Fig. 2] shows the mode frequencies as a function of Ma 
for ju — 2, q — 1, £^bv = ^-2 (same as for the case shown 
in Fig. 0, but now with a finite friction parameter: a = 0.2. 
The upper panel gives the mode structure in the same format 
as Fig. [2] The presence of friction removes the symmetry of 
the imaginary parts of a> (the growth rates) and also prohibits 
the mode merging at Ma — 0.55 at the onset of the Kelvin- 
Helmholtz instability for a = 0. As shown in detail in the lower 
panel of Fig. @] the first backward mode becomes unstable for 
Ma — 0.31 and its the growth rate shows a significant upturn 
for Ma £ 0.55, at which value the Kelvin-Helmholtz instability 
set in for a — 0. The second backward mode becomes unsta- 
ble at Ma - 1.45, but the system is by then already dominated 
by the first unstable mode. We note that, like in the unstrati- 
fied case, the friction-induced instability sets in earlier (i.e., for 
lower values of Ma) then the Kelvin-Helmholtz instability. 



4. Conclusions 

In the framework of the approximation of thin flux tubes, 
we have provided a unified linear treatment of the Kelvin- 
Helmholtz instability, the dissipative (friction-induced) insta- 
bility, and the undulatory (Parker-type) instability of a mag- 
netic flux tube with a field-aligned flow in a gravitationally 
stratified medium. Including dissipation effects by a Stokes- 
like friction term leads to overstability of transverse waves for 
flow velocities below the threshold of the Kelvin-Helmholtz in- 
stability. In the case of a stratified medium, the critical flow ve- 
locity can become arbitrarily low near the threshold for Parker 
instability. In other words, for a given flow velocity, there is al- 
ways a range of parameters (e.g., perturbation wavelength) that 
already give friction-induced instability before the Parker-type 
instability sets in. 

The existence of friction-induced instability (overstability 
of transverse waves) potentially has important consequences 
for magnetic flux storage in stellar convection zones. Toroidal 
magnetic flux tubes require faster internal rotation (equiva- 
lent to a longitudinal flow) in order to keep force equilibrium 
(Moreno-Insertis et al ., 1992). Although this flow is always 
sub-Alfvenic, we have seen that it still can lead to instability 
in a stratified medium. This will be studied in detail in the third 
paper of this series. 

To include the drag-induced dissipation in our linear analy- 
sis, we have replaced the aerodynamic drag term (quadratic in 
the transverse velocity) commonly used in treatments of thin 
flux tubes by a Stokes-like friction term. As we will show by 
numerical simulations in paper III, this does not affect the sta- 
bility criteria for the friction-induced instability derived here. 
In fact, the precise nature of the dissipation process is not cru- 
cial for the onset of the overst ability associated wi t h backward 



or neg ative-energy waves (cf., iJoarder et al. 
1998). 



1997; Tirrv et al 



Appendix A: Kelvin-Helmholtz instability for 
external flow 

The results of the stability analysis have to be independent of 
Galilean transformations of the frame of reference. Since, in 
our case, the effect of the external medium on the flux tube 
is described solely by the enhanced inertia for transverse ac- 
celeration (apart from the pressure balance condition), it is not 
completely straightforward to show that the result is Galilei in- 
variant. Therefore, we give the derivation in what follows. 

Assume a static straight equilibrium flux tube embedded in 
an external flow along the tube with u e = unfito- Consider a dis- 
placement, £„(*<)> Oi of a Langrangean mass element in the tube 
with normal acceleration £ n (so, t). The corresponding accelera- 
tion of a mass element outside the flux tube is given by 



d 
dr 



dr 



d7 [& + 



where 

d_ _ d_ _d_ 

dr dt dso 



(A.l) 



(A.2) 
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Dots indicate partial time derivatives and primes indicate par- 
tial derivatives with respect to arc length along the undisturbed 
tube. By using Eq. flA.lt and applying "actio = reactio", we find 
for the normal component of the linearized equation of motion 

Piofn - PiO«|^n = -PeO (in + 2t) £, + . (A. 3) 

Introducing fi — 1 + peu/PiO an d the ansatz £ n 00 exp(fco _ w^O 
into this equation, we obtain 

fi(b 2 + (p - 1) (^ 2 ^ - 2o M>) - /fc 2 ^ = . (A.4) 

We now carry out a Galilean transformation into the rest frame 
of the external medium. In this frame of reference, we have 
an internal longitudinal flow with velocity Moe t o = -fo^to- 
Inserting mo = -Vo and the transformation a> = ut-vok = cb+u^k 
into Eq. ( IA.4I ). we find 

/ioj 2 - 2u{)ka) + k 2 (iiy — UjA = , (A. 5) 

which is identical to Eq. ©. 

References 

Achterberg, A. 1982, A&A, 1 14, 233 

Achterberg, A. 1996a, A&A, 313, 1008 

Achterberg, A. 1996b, A&A, 313, 1016 

Andries, J. & Goossens, M. 2001a, A&A, 368, 1083 

Andries, J. & Goossens, M. 2001b, A&A, 375, 1 100 

Andries, J. & Goossens, M. 2002, Phys. Plasmas, 9, 2876 

Batchelor, G. K. 1967, An Introduction to Fluid Dynamics 

(Cambridge University Press) 
Cairns, R. A. 1979, J. Fluid Mech., 92, 1 
Chandrasekhar, S. 1961, Hydrodynamic and hydromagnetic 

stability (Oxford University Press, Oxford, UK) 
Cheng, J.-H. 1994, Chinese Astronomy and Astrophysics, 18, 

374 

D'Silva, S. Z. & Choudhuri, A. R. 1991, Sol. Phys., 136, 201 
Ferrari, A., Trussoni, E., & Zaninetti, L. 1981, MNRAS, 196, 
1051 

Ferriz-Mas, A. & Schiissler, M. 1993, Geophys. Astrophys. 

Fluid Dyn., 72, 209 
Ferriz-Mas, A. & Schiissler, M. 1995, Geophys. Astrophys. 

Fluid Dyn., 81,233 
Goossens, M., Hollweg, J. V., & Sakurai, T. 1992, Sol. Phys., 

138,233 

Hollweg, J. V., Yang, G, Cadez, V. M., & Gakovic, B. 1990, 
ApJ, 349, 335 

Holzwarth, V. 2002, PhD thesis, University of Gottingen, 

Germany, http://webdoc.sub.gwdg.de/diss/2002/holzwarth 
Joarder, P. S., Nakariakov, V. M, & Roberts, B. 1997, 

Sol. Phys., 176, 285 
Kolesnikov, F, Biinte, M., Schmitt, D., & Schiissler, M. 2004, 

A&A, 420, 737 
Lashmore-Davies, C. N. 2005, J. Plasma Phys., 71, 101 
Moreno-Insertis, F, Schiissler, M., & Ferriz-Mas, A. 1992, 

A&A, 264, 686 

Moreno-Insertis, F, Schiissler, M., & Ferriz-Mas, A. 1996, 

A&A, 312, 317 
Nakariakov, V. M. & Roberts, B. 1995, Sol. Phys., 159, 213 



Osin, A., Volin, S., & Ulmschneider, P. 1999, A&A, 351, 359 

Parker, E. N. 1964, ApJ, 139, 690 

Parker, E. N. 1975, ApJ, 198, 205 

Rae, I. C. 1983, A&A, 126, 209 

Ray, T. P. 1981, MNRAS, 196, 195 

Ryutova, M. P. 1988, J. Exp. Theor. Phys./JETP, 67, 1594 

Ryutova, M. P. 1990, in Solar Photosphere: Structure, 

Convection and Magnetic Fields, IAU Symp. 138, ed. J. O. 

Stenflo (Dordrecht: Kluwer Academic Publishers), 229-249 
Satya Narayanan, A. & Somasundaram, K. 1982, Ap&SS, 84, 

247 

Schmitt, D. 1998, Geophys. Astrophys. Fluid Dyn., 89, 75 
Schiissler, M. 1990, Habilitationsschrift, Universitat 

Gottingen, ArXiv Astrophysics e-prints: astro-ph/0506050 
Schiissler, M. & Ferriz-Mas, A. 2007, A&A, 463, 23 
Spruit, H. C. 1981, A&A, 102, 129 

Spruit, H. C. & van Ballegooijen, A. A. 1982, A&A, 106, 58 
Tiny, W. J., Cadez, V. M., Erdelyi, R., & Goossens, M. 1998, 

A&A, 332, 786 
van Ballegooijen, A. A. 1983, A&A, 275 
Walker, A. D. M. 2000, J. Plasma Phys., 63, 203 



